Identifying wave packet fractional revivals by means of information entropy 
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Wave packet fractional revivals is a relevant feature in the long time scale evolution of a wide 
range of physical systems, including atoms, molecules and nonlinear systems. We show that the 
sum of information entropies in both position and momentum conjugate spaces is an indicator of 
fractional revivals by analyzing three different model systems: (i) the infinite square well, (ii) a 
particle bouncing vertically against a wall in a gravitational field, and (in) the vibrational dynamics 
of hydrogen iodide molecules. This description in terms of information entropies complements the 
usual one in terms of the autocorrelation function. 
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The phenomenon of quantum wave packet revivals has 
received wide attention over the last years. It has been 
investigated theoretically in atomic and molecular quan- 
tum systems and observed experimentally in, among 
others, Rydberg wave packets in atoms and molecules, 
molecular vibrational states, and Bose-Einstein conden- 
sates 3- Revivals occur when a wave packet solution 
of the Schrodinger equation evolves in time to a state 
that closely reproduces its initial waveform. Fractional 
revivals appear as the temporal formation of structures 
that are given by a superposition of shifted and rephased 
initial wave packets [3|, y, |5|. It has been shown 
that the relevant time scales of wave function evolution 
are contained in the coefficients of the Taylor series of the 
energy spectrum. En, around the energy Eng correspond- 
ing to the peak of the initial wave packet. More precisely, 
the second-, third-, and fourth-order terms in this expan- 
sion are associated with, respectively, the classical period 
of motion T^i, the quantum revival time-scale Tret,, and 
the so-called superrevival time. Fractional revival times 
can be given in terms of the quantum revival time-scale 
by i = pTrev/q, with p and q mutually prime 0], and 
are usually analyzed using the autocorrelation function 
A{t), which is the overlap between the initial and the 
time-evolving wave packet [Tj. At certain fractional re- 
vivals, however, the autocorrelation function may be of 
limited help since the wave packet reforms itself, possibly 
into a scaled copy of its original shape, in a location that 
does not generally coincide with its initial position. An 
expectation value analysis of wave packet evolution has 
been recently proposed by some authors^but it misses to 
fully detect the fractional revivals jH, HI, Q . 

In this Letter we study the wave packet dynamics by 
means of the sum of the information entropies of the 
probability density of the wave packet, in both position 
and momentum spaces. We shall show that it provides 
a natural framework for fractional revival phenomena. 
The position-space information entropy measures the un- 
certainty in the localization of the particle in space, so 
the lower is this entropy the more concentrated is the 



wave function, the smaller is the uncertainty, and the 
higher is the accuracy in predicting the localization of 
the particle. Momentum-space entropy measures the un- 
certainty in predicting the momentum of the particle. 
Thus, information entropy gives an account of the spread- 
ing (high entropy values) and the regenerating (low en- 
tropy values) of initially well localized wave packets dur- 
ing the time evolution. Moreover, if p{x) = \ip{x)\'^ 
and 7(p) — |0(p)P are respectively the probability den- 
sities in position and momentum spaces (where ip and 
(j) are the position and momentum wave packets), the 
uncertainty relation for the information entropy implies 
5p -I- 5^ > 1 -I- Invr, where Sp = — J p{x) hip{x)dx and, 
analogously, = — J 7(p) hij{p)dp. This inequality is a 
generalization of the standard variance-based Heisenberg 
uncertainty relation 10] . It is satisfied as a strict equality 
only for Gaussian wave packets and bounds from below 
the sum of the entropies to l-|-ln7r. During the evolution 
of a Gaussian wave packet, the entropy sum decreases at 
the revival times to reach the above value, which plays 
a similar role to unity in the autocorrelation function. 
Furthermore, the formation of a number of 'minipackets' 
of the original packet, i.e. fractional revivals of the wave 
function, will correspond to the relative minima of the 
total entropy. Notice that it is the sum of the entropies 
that is employed as an indicator of the fractional revivals, 
and not either of them separately, because only the sum 
embraces both the configurational and the motion as- 
pects of the wave packet dynamics. In this context, we 
point out that the sum of entropies for the phase and 
photon number has been used to study the formation of 
macroscopic quantum superposition states from an ini- 
tially coherent state via interaction with a Kerr medium 
[ll| . Additionally, we note that a finite difference eigen- 
value method has been recently derived from which the 
various orders of revivals can be directly calculated rather 
than searching for them [12:]. 

In what follows we investigate the time evolution of 
wave packets in three one-dimensional, model systems 
that exhibit fractional revival behavior, namely, the fa- 
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miliar infinite square-well, the quantum 'bouncer', that 
is, a quantum particle bouncing on a hard surface under 
the influence of gravity, and a superposition of molecular 
wave packets in a Morse potential describing the vibra- 
tions of hydrogen iodide molecules. The infinite square- 
well is especially well suited for the understanding of frac- 
tional revival phenomena because of its amenability to 
analytical treatment, and has been analyzed by several 
authors (see, for instance, 0, H, [H ) . On the other hand, 
the quantum bouncer is the quantum version of a fa- 
miliar classical system and has also been studied in the 
context of wave packet propagation and revival phenom- 
ena Recently, gravitational quantum bouncers 
have been realized using neutrons and atomic clouds 
[lit I hence endowing this system with great physical sig- 
nificance. Finally, molecular wave packets provide a real- 
istic scenario for the assessment of the entropy approach. 
Revivals and fractional revivals have been observed ex- 
perimentally in wave packets involving vibrational levels 
and can beprobed by random-phase fluorescence inter- 
ferometry [2[l9l| 

Consider an infinite potential- well defined as V{x) =0 
for < a; < L and V{x) = +oo otherwise. The time- 
dependent wave function for a localized quantum wave 
packet is expanded as a one-dimensional superposition 
of energy eigenstates as 
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where Un{x) represent the normalized eigenstates 
and En the corresponding eigenvalues, Un{x) — 
y/2/Lsm (n-Kx/L) , En — n^/i^7r^/2mi^. Following the 
customary procedure, the classical period and the re- 
vival time can be computed as Td = 2mLp'/fmn and 
Trev = 'imL'^ /hn, respectively. As an example, it is easy 
to see by direct substitution in ^ that ip{L—x, Trev/"^) = 
—•ip{x, 0), so at time t = Trev/'2 a copy of the initial state 
reforms itself, reflected around the center of the well 4]. 

We shall consider an initial Gaussian wave packet with 
a width CT, centered at a position xq and with a mo- 
mentum pq, ip{x,0) — exp[— (x — xo)^ /2a^h^ + ipa{x — 
xq ) / fi.] / ■\/ ahy/iT. Assuming that the integration region 
can be extended to the whole real axis, the expan- 
sion coefficients can be approximated with high accuracy 
by an analytic expression Q. To calculate the corre- 
sponding time dependent, momentum wave function we 
use the Fourier transform of the equation ([T|), and the 
momentum-space normalized eigenstates 



7tL 



9 

Pn 



i-iy 



jpL/h 



1 



(2) 



Without loss of generality, we shall henceforth take 2m — 
h — L = 1, and a = 1/10 for the initial wave packet. 

In Fig. [1] the sum of entropies, Sp[t) + S-y{t), and 
the autocorrelation function, are shown for an ini- 

tial wave packet with xq = 0.5 and po = 4007r. At early 




FIG. 1: Time dependence of \A{t)f and Sp{t) + S-y{t) for 
an initial Gaussian wave packet with xo = L/2, po = 4007r, 
and cr = 1/10 in an infinite square- well. A) Long-time depen- 
dence. The main fractional revivals are indicated by vertical 
dotted-lines. B) First classical periods of motion. Quantum 
wave packets are represented by solid lines and their classical 
component by dashed-lines. 



times, the Gaussian wave packet evolves quasi-classically, 
but in a few periods the quantum and classical wave 
packet trajectories start moving apart (Fig. [Tj3), the 
classical component of the wave function being defined 
as ipci{x, t) = a„M„(a;)e~*^'^"*/^'=' For longer time 
scales, a large amplitude modulation is superimposed on 
the quasiperiodic oscillations (Fig. [H'V). In this long-time 
regime, the wave packet initially spreads and delocalizes 
while undergoing a sequence of fractional revivals with 
the creation of sub-packets, each of them similar to the 
initial one so that the sum of the entropies reaches a 
relative minimum and the modulus of the autocorrela- 
tion function a relative maximum. The most important 
fractional revivals are denoted by vertical dashed-lines. 
It can be observed that the identification of some frac- 
tional revivals from the autocorrelation function is not 
so clear-cut as compared with the entropy approach (see, 
for example, the cases t = pTrev/l with q = 7 or 9). 

We have also investigated the interesting initial con- 
dition xo = 0.8L and po = 0, for which there is no 
classical periodic motion. Snapshots of the numerical 
simulation of the position-space probability density are 
given in Fig. [5] at several times. It is apparent from 
the bottom panel of Fig. [3] that the sum of entropies 
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FIG. 2: Snapshots of the probabihty density in position-space 
for an initial Gaussian wave packet with xq = 0.8L, po — 
0, and (J = 1/10 in an infinite square- well and at different 
fractional revival times. 



FIG. 4: Time dependence of (top panel) and (bottom 

panel) Sp{t) + S~,{t) for an initial Gaussian wave packet with 
zo — 100, po = 0, and ct = 1 in a quantum bouncer. The main 
fractional revivals are indicated by vertical dashed-lines. 
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FIG. 3: Time dependence of (top panel) |A(f)|^ and (bottom 
panel) Sp{t) + S-y{t) for an initial Gaussian wave packet in an 
infinite square well. Parameters as in Fig. [S] 
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FIG. 5: Time dependence of (top panel) |A(f)|^ and (bottom 
panel) Sp{t) + S^{t) for an initial superposition of molecular 
wave packet with Gaussian weights in a Morse potential. The 
main fractional revivals are indicated by vertical dashed-lines. 



has a minimum at the main fractional revivals, denoted 
by vertical dashed-lines, except for the case t = Trev/7 
where the position-space probability density has a shape 
compatible with a collapsed wave function (Fig. [2|). The 
autocorrelation function, as plotted in the top panel of 
Fig. [31 fails to show the fractional revivals occurring 
at, for example, t/Trev = 1/6, 1/8, and 3/10. This dis- 
crepancy can be traced back to the fact that information 
entropies take into account individual minigaussian pack- 
ets independently of their relative positions, whereas the 
autocorrelation function depends on the relative position 
between the initial wave packet and the evolved one. 

Next, we study the behavior of these same quantities 
in other physical situations. Consider a particle bounc- 
ing on a hard surface under the influence of gravity, that 



is, a particle in a potential V{z) — mgz, if z > and 
V{z) ^ +00 otherwise. Upon introducing the character- 

1/3 

istic gravitational length Ig = (h/2gm^) and defining 
z' — z/lg and E' = E/mglg, as in [20|, the eigenfunc- 
tions and eigenvalues are given by E!^ = ZmUn{z') = 
AfnAi{z' — Zn) with n = 1,2,3, .. ., where Ai{z) is the 
Airy function, ~Zn denotes its zeros, and Mn is the Un{z') 
normalization factor. Very accurate analytic approxima- 
tions to Zn and Mn can be found in [2(]| . Consider now 
an initial Gaussian wave packet localized at a height zq 
above the floor, with a width a and an initial momen- 
tum Po — 0. The corresponding coefficients in Eq. ^ 
can be obtained analytically [l^l , and the classical period 
and the revival time are Td = and T^ev = 4zQ/7r, 

respectively. 
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The temporal evolution of the total entropy and of 
the autocorrelation function was computed numerically 
for the initial conditions Zq = 100, a ~ 1, and Po ~ 
(see Fig. For that, the corresponding wave packet in 
momentum-space was obtained numerically by the fast 
Fourier transform method. One can see in the bottom 
panel of Fig. [4] that the reformation of subpackets at 
fractional revivals is captured by the successive relative 
minima of Sp{t) + S-y{t). A similar information is pro- 
vided by the, somewhat less clear, sequence of relative 
maxima of (top panel of Fig. 2]). 

Last, we address the case of a diatomic molecule, the 
vibrational dynamics of which is known to be well ap- 
proximated by the Morse potential, V{x) = D{e~'^^^ — 
2e~^^). Here, x — r/vQ — 1, r is the internuclear dis- 
tance, ro is the equilibrium bond distance, D is the dis- 
sociation energy and /3 is a range parameter. Defining 
A = \/2jiDro/Pfi and s = 2\\J —E j D , the associated 
eigenvalues and eigenfunctions for bounded states can 
be written as [2l| ^'^(O = ^ En = 
-D{X - n - l/2)VA^ with ^ = 2Ae"'3^, < < c», 
and n = 0, 1, [A — 1/2], [x] being the integer part of x. 
^n(C) ^re the Laguerre polynomials and is a normal- 
ization factor (see [2l|). 

We shall consider the HI molecule, for which /3 = 
2.07932, D = 0.1125 a.u., tq = 3.04159 a.u., with the 
reduced mass /i = 1819.99 a.u. The number of bound 
states is [A — 1/2] + 1 = 30, and the classical period 
and the revival time are given by Td = Trev/i'2X — 1) 
and Trev = 2ttX^/D, respectively. Notice that in this 
case the autocorrelation function turns out to be sym- 
metric about Trei,/2 [2^ . As an initial wave packet, we 
shall take a superposition of Morse eigenstates assum- 
ing a gaussian population of vibrational levels |c„p = 
exp ((n — no)^/(T) /(ttct)^/^, with no = 7, a = 3. As in 
the quantum bouncer case, the fast Fourier method is 
used to compute the entropy in the momentum space. 
These system parameters lead to the entropy and auto- 
correlation function depicted in Fig. [5l where the oc- 
currance of fractional revivals at t = 1/6,1/3,2/3, and 
5/6 is transparent at first sight in the bottom panel, in 
contrast with the information provided by the autocor- 
relation function (top panel). 

In conclusion, we have found that the manifestation 
of fractional revivals in the long time-scale evolution of 
quantum wave packets reflects in the sum of information 
entropies in conjugate spaces, which clearly shows rela- 
tive minima at the fractional revival times, thus provid- 
ing a useful tool to reckon with for visualizing fractional 
revivals, complementary to the conventional autocorre- 
lation function. These minima appear independently of 
the relative position of the subpackets that configure the 
wave packet at the fractional revival time, and are rigor- 
ously bounded from below. It would be very interesting 
to carry out a more systematic study extending this ap- 



proach to systems with two or more quantum numbers. 
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